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COMPARISON OF FLUTTER CALCULATIONS USING 
VARIOUS AERODYNAMIC COEFFICIENTS WITH EXPERIMENTAL 
RESULTS FOR SOME RECTANGULAR CANTILEVER 
WINGS AT MACH NUMBER 1.3 


By Herbert C. Nelson and Ruby A. Rainey 
SUMMARY 


A general Rayleigh analysis is used as a basis for developing four 
methods of flutter analysis that are applied to twelve low-aspect-ratio 
wings. These wings were previously tested at a Mach number of 1.3 by 
progressively varying certain wing parameters until flutter occurred. 
They were rectangular in plan form and had aspect ratios between 3.00 
and 4.55. The four methods of flutter analysis used are: section coef- 
ficients for harmonically pitching and translating rectangular wings in 
& Rayleigh type of analysis, two-dimensional coefficients in a Rayleigh 
type of analysis, total coefficients for harmonically pitching and trans- 
lating rectangular wings in a representative-section analysis, and two- 
dimensional coefficients in a representative-section analysis. Each of 
the four methods involved two degrees of freedom, namely, first bending 
and first torsion of a cantilever wing. 


The analytical results are compared with the previously obtained 
experimental values. The comparison indicates that the use of section 
aerodynamic coefficients derived on the basis of three-dimensional fiow 
leads to & significant improvement in the correlation of theory and 
experiment. 


INTRODUCTION 


The problem of theoretically determining the flutter characteristics 
of unswept wings of low aspect ratio in supersonic Plow has become of 
increased interest. Most of the previous analytical work on this problem 
has been based on air-force and moment coefficients for two-dimensional 
supersonic flow, such as those tabulated in reference 1. For example, 
reference 2 presents the results obtained at a Mach number of 1.5, by 
using two-dimensional coefficients in a representative-section type of 
flutter analysis, for twelve unswept wings with aspect ratios ranging 
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from 2.00 to 4.55. As explained in reference 2, these wings were also 
tested at a Mach number of 1.3 by progressively shifting their centers 
of gravity and elastic axes and modifying their bending and torsional 
frequencies until flutter occurred. A comparison of the calculated and 
experimental results showed that in the majority of cases the calculated 
flutter speeds were considerably below the experimental flutter speed. 
This discrepancy suggests in part that, at least in the low supersonic 
speed range, two-dimensional coefficients are inadequate and more real- 
istic aerodynamic coefficients should be used in the flutter analysis 

of unswept low-aspect-ratio wings. 


In reference 5, streamwise section and total air-force and moment 
coefficients expanded to the seventh power of the frequency of oscilla- 
tion were developed for harmonically pitching and translating rectangu- 
lar wings moving at supersonic speed. The section coefficients were 
used in a Rayleigh type of flutter analysis to calculate the flutter 
speeds of a rectangular wing of aspect ratio 4.53 at several Mach num- 
bers in the low supersonic speed range. For comparison the wing was 
also analyzed by using the two-dimensional coefficients of reference 1 
in a Rayleigh type of analysis. Examination of the results showed the 
flutter speeds based on the rectangular-wing section coefficients to be 
higher than those based on two-dimensional coefficients, particularly 
at the lower Mach numbers. Application of a Rayleigh type of analysis 
involving the section coefficients of reference ^ to the wings of ref- 
erence 2 might therefore be expected to yield a better correlation between 
theory and experiment than was obtained in reference 2. 


Also of interest is reference } in which a comparison is made between 
flutter results obtained by using two-dimensional coefficients in a 
representative-section type of analysis and total coefficients for rec- 
tangular wings in the same type of analysis. For wing parameters in the 
range of those given in reference 2, reference lh also shows an increase 
in calculated flutter speed resulting from the use of finite-wing 
coefficients. | 


Tn the present paper four methods of analysis are applied to the 
twelve wings of reference 2 and the results are compared with the exper- 
imental results in reference 2. These four methods of flutter analysis 
are: section coefficients for а pitching and translating rectangular wing 
in а Rayleigh type of analysis, two-dimensional coefficients in а Rayleigh 
type of analysis, total coefficients for a pitching and translating rec- 
tangular wing in a representative-section type of analysis, and two- 
dimensional coefficients in a representative-section type of analysis. 
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SYMBOLS 


aspect ratio, s/b 
one-half chord 
speed of sound in undisturbed medium 


first bending and first torsion damping coefficients, 
respectively (see ch. IX of ref. 5) 


vertical displacement of axis of rotation Χρ, positive 
downward. 


generalized coordinate in bending degree of freedom, hoe iot 
bending amplitude at tip of wing 
reduced frequency, ab/V 


coefficients of section lift and moment, respectively, 
associated with mode shape Zh 


coefficients of section lift and moment, respectively, 
associated with mode shape Zæ 


components of section force and moment coefficients, respec- 
tively, for rectangular wing (see ref: 4) in equation uU 
and for two-dimensional wing (see ref. 1) in equation (8); 
1 = 1, 2, 5, and k 


components of total force and moment coefficients, respec- 
tively, for rectangular wing (see ref. 5); i = 1, 2, 3, 
and Н 


Mach number, V/c 


aerodynamic section moment on wing about axis of rotation Xp, 
positive leading edge up 


aerodynamic section normal force, positive downward 


nondimensional radius of gyration of wing section about 


elastic axis, үт /tib= where I, is mass moment of 


inertia per unit span about elastic axis and m is mass 
of wing per unit span 


td 
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one-half span of wing 
time 
velocity of flow 


nondimensional chordwise coordinate measured from leading 
edge, referred to wing chord Zb 


location of center of gravity of wing measured from elastic 
axis (see ref. 1) 


chordwise position of axis of rotation of wing (elastic axis) 


nondimensional spanwise coordinate measured from midspan of 
wing, referred to wing half-span в 


first bending mode shape of wing 

first torsion mode shape of wing 

angle of attack, positive leading edge up 

Layt 


generalized coordinate in torsion degree of freedom, age 


torsion amplitude at tip of wing 


density parameter, xpb?/m 
density in undisturbed medium 
frequency of oscillation at flutter 


first bending frequency of wing 


first torsion frequency of wing 


METHODS OF FLUTTER ANALYSTS 


Rayleigh Analysis 


General considerations.- The wings to be analyzed are rectangular 


in plan form and were tested as cantilevers in the Langley supersonic 


NACA TN 3301 2 


flutter apparatus (a 9-inch by 18-inch supersonic drawdown tunnel). In 
& Rayleigh type of analysis of such wings, the bending component of the 
flutter mode can be approximated by the first bending mode of a uniform 
cantilever wing and the torsion component by the first torsion mode. 
The flutter determinant is then formed and is solved for the flutter 
condition. (A detailed discussion of the Rayleigh type of analysis as 
applied to flutter may be found in ch. IX of ref. 5.) 


The bending component h and the torsion component a of the 
flutter mode may be written as 


h(y,t) = ZM (y)h(t) 


(1) 


a(y,t) = 9(у)а(%) 


where y is the nondimensional coordinate shown in figure 1, Z and 
Z, are the first bending and first torsion mode shapes shown in fig- 


ure 2, and ћ and G are the generalized coordinates in the bending 
and torsion degrees of freedom, respectively. The section aerodynamic 
force or aerodynamic force per unit span, positive downward, associated 
with equations (1) may be written as 


P = ъа? [n GR + 1а(у)&) (2) 


and the section moment, positive leading edge up, about the arbitrary 
axis of rotation x = Χο may be written as 


Μα = gts [αρ (y) + ma(y)a| (5) 


where œ is the frequency of oscillation, Ὁ is the one-half chord о? 
the wing, I. and m are complex coefficients of the lift and moment 


associated with the mode Zp, and lą and m, are the complex coeffi- 
cients of the lift and moment associated with the mode Ζα. Each of the 
aerodynamic coefficients lh; Mh; la, and mg, in addition to being а 


funetion of the spanwise variable y, is & function of Mach number M 
&nd reduced frequency К = bo/V. Although these coefficients may be 
taken to apply at either subsonic or supersonic speed, the present paper 
is concerned only with the supersonic speed range. 
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The equilibrium equations at flutter may be obtained by setting up 
the potential and kinetic energies and the work of the applied forces, 
introducing equations (1), (2), and (3) and the mass and stiffness prop- 
erties of the wing, and then applying Lagrange's dynamical equation, as 
shown in chapter IX of reference 5. From the equilibrium equations & 
flutter determinant may be obtained in the form 


Ahh Ang, 
= 0 (4) 
Aah Aag, 


where the determinant elements are 


BOE Ίμζῃ dy 


Ahh = 
h T «ο 
tw =a Í адабу - κ É 102 dy 
(5) 
СЕК 
Ав = x | Palm ay - Z x J mhg, dy 


2 1 
к= seb - (P| P =s а [ пода, ду 


From equations (4) and (5) four methods of analysis are obtained by using 
various approximations in evaluating the integrals of equation (5). 


Section coefficients for rectangular wing.- The following approxi- 
mate expressions for the section coefficients are employed: 
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th Zy (Ia + ilo) 


lg, = Za [E + 11) 


(6) 


mh Zr, (M; + мо) 


My = Z, (M5 + м) 


where Τη and M; (1 = 1, 2, 3, and 4) are the components of the sec- 


tion coefficients given in reference 5 for & rectangular wing in super- 
sonic flow oscillating harmonically as a rigid body in pitch and verti- 
cal translation. (A preliminary unpublished analysis, based on parabolic 
bending of a rectangular wing which closely resembles the mode shape Zh, 


suggests that the results obtained by using the distributions of lift and 
moment for the mode shapes Z, and Z would be nearly identical to the 


results obtained by using the approximate distributions given by eqs. (6), 
when multiplied by the mode shape Zh or Z, and integrated in the man- 


ner required in eqs. (5).) Upon substituting equations (6) into equa- 
tions (5), the determinant elements of equation (4) become 


mn = b -(2) | EFIE 
мш = x0 [аа-а f 


(25 + 1n) Zi Ze, dy 


(7) 
1 
ка = |” Zaza ἂν - $ e J" (щ + ne) as 


кш, = zala - (82) Г ау [P (5 + mayas 
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The uncoupled first bending mode shape ση and the first torsion mode 
shape 9, needed for the evaluation of the integrals of equations (T) 


are shown in figure 2. The integrals of equations (7) containing only 
mode shapes can be evaluated to give 


Г 7.28у = 0.25 


1 
I 2124, dy = 0.5357 


A numerical method for evaluating the integrals of equation (7) involving 
the aerodynamic coefficients L; and M, is given in appendix B of ref- 
erence 5. (In using ref. 2, note that the spanwise coordinate у о? the 
present paper and the spanwise coordinate € of the reference paper are 
related by у = 1 - t.) 


Coefficients for two-dimensional wing.- If two-dimensional air-force 
&nd moment coefficients &re used in place of the section coefficients of 
reference 5, the force апа moment coefficients in equations (T) appear 
as constants in the integrals and can be factored from under the integral 
signs, and the determinant elements of equation (4) become 


wi P= (BY Eee ena] at 
ha [κω - 8 a(t + atu) | Ζοζα ἂν 
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where 1} and Μι (1 = 1, 2, 5, and h) now refer to components of 
be two-dimensional coefficients, such as those tabulated in reference 1. 


Representative-Section Analysis 


Total coefficients for rec ar wing.- By applying mean-value ` 
theory to the integrals in equation © and, in the process, by assuming 
the representative section to be the same for all integrals involved, 
the determinant elements can be written as 


| 

HP 

l 
— 
E 

1 
Aff 


k Г (ta + по) ву (2,2), 
eee (a+ gen. 


where the subscript r denotes evaluation at a representative spanwise 
station у = г. Since the quantities having the subscript г cancel 
in the solution of equation (N), equations (9) may be rewritten ав 


ha = 1 - (8) - $ E + οἷς 


Ang = Χα - + (is + ify) (10) 


Ag = xo = (Rh + 382) 


. Ao, = xh - СЙ τ. + 2) 
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l 
where L, -Γ Ly dy and Μι z: M; dy (1-21, 2, 3, and h) are the 


components of the total force апа moment coefficients for а pitching and 
translating rectangular wing given in reference 5. 


Coefficients for two-dimensional wing.- If infinite aspect ratio is 
substituted into the aerodynamic coefficients of equation (10) (see 
ref. 5), the determinant elements can be written as 


(11) 


where, as in equation (8), I4 and M, refer to components of two- 
dimensional coefficients, such as those tabulated in reference l. 


Solution of Flutter Determinant 


The flutter condition is determined from the nontrivial solution 
of equation (4) obtained by using as determinant elements the various 
approximate forms of equations (5) given by equations (7), (8), (10), 
and (11). This condition, which requires that the real and imaginary 
parts of equation (4) vanish simultaneously for the same set of aero- 
dynamic and wing parameters, may be obtained by various means (see 
ch. XIII of ref. 5). 


In the present paper the ratio (y/o in equation (4) is replaced 
by the equivalent quantity (ох, fa) (ax [oo ) Then, for a particular wing 
and Mach number, for which values of M, к, Χο, Χα, rg2, and ах, Jo, 
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are specified, equation (4) contains the two unknown parameters 0/0 


апа k = bov. The reduced frequency к (upon which the various aero- 
dynamic terms are dependent) is varied until the same value of ay/w 


is obtained from both the real and imaginary parts of equation (4). ‘This 
is the required condition and yields the values of k and а/о at 


flutter and consequently the flutter-speed coefficients V/ bax, for the 
wing at the selected value of M. 


APPLICATION AND DISCUSSION OF RESULTS 


The four methods of analysis outlined in the previous section were 
applied to the twelve wings of reference 2. The wing parameters needed 
in these analyses and a description of each wing profile, obtained from 
reference 2, are given in table I of the present paper. Тһе flutter 
parameters, v/ bo and w/t, and consequently v/ бам, calculated by 


these methods are listed in table II. For comparison table II also 
includes the experimentally determined flutter parameters given in 
reference 2. 


In figure 5 the data of table II are plotted in line-graph form. 
The line-graph method of plotting is employed to achieve a separation 
of the data and ease of comparison not otherwise obtained because of 
the insufficient range of variation of the different wing parameters. 
Also shown in figure 5, as flagged points, are the analytical results 
of reference 2. These results were obtained by the last method of the 
previous section (two-dimensional coefficients in a representative- 
section analysis) but included structural damping. Structural damping 
could also have been included in the calculations of the present paper by 


replacing ax? by an? (1 + ign) and. а2 by ay? (1 t iga)» where өр 
is the damping coefficient in bending and gx ів the damping coefficient 


in torsion, in the methods discussed previously. Since damping was not 
included, the calculations of reference 2 may serve to indicate the 
effect the inclusion of damping would have on the calculations of the 
present paper. 


Figure 5(а) shows a comparison for each wing of the values of 
reduced flutter speed V/bæ (reciprocal of reduced frequency k), 
obtained by the four methods of analysis. The results obtained by 
using finite-wing section coefficients in a Rayleigh analysis are 
closest to experiment in all twelve of the cases treated. 


In figure 3(b) values of the ratio of flutter frequency to tor- 
sional frequency w/c, are compared for the various wing models. As 
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may be noted in the figure, for five of the wing models (A-1, B-1, C-1, 
C-2, and F-1) comparatively large differences between the theoretical 
and experimental values of оја, exist; these differences would probably 


be reduced by the inclusion of more degrees of freedom in the various 
analyses. However, it may be seen by comparing the flagged and unflsgged 
right triangles in figure 3(b) that the inclusion of structural damping 
in the Rayleigh analysis involving section coefficients for a rectangu- 
lar wing may sufficiently reduce the differences between experiment and 
theory. 


Figure 5(c) shows a comparison of the values of flutter-speed coef- 
ficient V/ bax, determined for the various wing models from the data 


presented in figures 3(a) and 3(b). Тһе results of using the Rayleigh 
analysis involving section coefficients for rectangular wings &re 
closest to experiment in the majority of the cases treated, that is, 
except for models A-1, C-1, and C-2. The section-coefficient results 
in these cases are above the experimental values (nonconservative). 
Inclusion of more modes in the analysis may relieve this situation. 


Also of interest in the present comparison are the curves of V/ dary, 


calculated in reference 5 for model B-1 of table I in the Mach number 
range 10/9 < M < 10/6 by the first two methods of the previous section, 


that is, section coefficients for a nondeformable rectangular wing in a 
Rayleigh type of analysis and two-dimensional coefficients in the same 
type of analysis. These curves, taken from figure 12 of reference 3, 
are shown in figure h. The main feature of these curves, as pointed 

out in the reference paper, is that the use of finite-wing coefficients 
is very influential at Mach numbers near unity but, as would be expected, 
becomes less so as the Mach number is increased. At М = 10/9, for the 
particular wing analyzed, the flutter speed obtained by using two- 
dimensional coefficients is about 62 percent of that obtained by using 
rectangular-wing section coefficients, whereas at М = 10/6 it is about 
95 percent. For comparison at М = 1.5 the experimental value for 
model B-1 and the results of using two-dimensional coefficients with 
and without structural damping and total rectangular-wing coefficients 
in а representative-section analysis are included in figure h. Тһе 
values plotted in figure 4 at М = 1.5 are, of course, also given for 
model B-1 in figure 5(с). Ав may be noted in figure h, the result 
obtained by using rectangular-wing section coefficients in a Rayleigh 
analysis is in excellent agreement with experiment. 


CONCLUDING REMARKS 


The results of applying four methods of flutter analysis to a 
series of twelve wings have been presented and discussed. The wings in 
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question, which were fluttered previously &t а Mach number of 1.5 in the 
Langley supersonic flutter apparatus, had aspect ratios ranging from 3.00 
to 4.55 and various profile shapes, masses, and stiffness properties. 

The four methods of analysis, which are derivable from a general Rayleigh 
type of analysis, are: section coefficients for a pitching and trans- 
lating wing in a Rayleigh type of analysis, two-dimensional coefficients 
in a Rayleigh type of analysis, total coefficients for a pitching and 
translating rectangular wing in a representative-section analysis, and 
two-dimensional coefficients in a representative-section analysis. Each 
of the four analyses involved two degrees of freedom, namely, first 
bending and first torsion of a cantilever wing. The section and total 
aerodynamic coefficients for rectangular wings that were used are those 
that were developed, for wing pitching and vertical translation, to the 
seventh power of the frequency in NACA TN 4076. 


Previous analyses of the flutter of unswept wings of low aspect 
ratio in supersonic flow have customarily involved the use of aerodynamic 
coefficients for two-dimensional flow. The present paper shows that the 
use of aerodynamic coefficients for rigid-body motions of a wing, namely 
pitching and vertical translation, derived on the basis of three- 
dimensional flow leads, at least in the low supersonic speed range, to 
& significant improvement in the correlation of theory and experiment. 


Langley Aeronautical Laboratory, 
Natjonal Advisory Committee for Aeronautics, 
Langley Field, Va., August 13, 1954. 
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TABLE II == 
COMPARISON OF CALCULATED AND EXPERIMENTAL FLUTTER PARAMETERS 
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Wing root 


| Midchord 
Leading ΠΤ 5 Trailing 
edge 2 edge 
Axis of X 


rotation 
ΧΞχο 
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(b) Section A-A. 


Figure l.- Illustration of coordinate system and two degrees of 
freedom a and n. 
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Figure 2.- Uncoupled first bending mode shape Zp and first torsion 
mode shape Z for a uniform cantilever wing. 
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Figure 5.- Comparison of calculated and experimental flutter character- w 
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istics for twelve wing models described in table I. 
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(b) Ratio of flutter frequency to torsional frequency, ur. 
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Figure 5.- Continued. 
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——— Rect. wing section coeff.l Rayleigh 
-— — — Two-dimensional coeff. onalysis 


A Rect. wing total coeff. Representative 
4 Two-dimensional coeff.| section analysis 


4 Calc. results of ref. 2, 9, and q, #0 


Figure 4.- Flutter-speed coefficients plotted against Mach number for 
model B-l of table I. (A = 4.533 £ = 95.5; xo = 0.342; χα = 0.250; 
Re = 0.59; and = = 0.585.) 
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